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ABSTRACT 

In this paper we perform a new estimate of the high energy neutrinos expected from GRBs associated 
with the first generation of stars in hght of new models and constraints on the epoch of rcionization and 
a more detailed evaluation of the neutrino emission yields. We also compare the diffuse high energy 
neutrino background from Population HI stars with the one from "ordinary stars" (Population II) , as 
estimated consistently within the same cosmological and astrophysical assumptions. In disagreement 
with previous literature, we find that high energy neutrinos from Population HI stars will not be 
observable at current or near- future neutrino telescopes, falling below both the sensitivity of a km'^ 
telescope and the atmospheric neutrino background, also under the most optimistic predictions for 
the GRB rate. This rules them out as a viable diagnostic tool for these still elusive metal-free stars. 
Subject headings: stars: early-type - gamma rays: bursts, neutrinos - cosmology: theory 



1. INTRODUCTION 

The first generation of stars (Population HI or PopIII 
stars) born after the collapse of the very first structures 
in the universe is puzzling the astrophysical community 
since long time. The vanishing metall icity of the uni - 
verse at the epoch of their formation (jlocco et al. '07l ) 
is supposed to give them peculiar properties, all aris- 
ing from their very high char acteristic mass d ue to the 
peculiar cooling of the cloud (jAbel et al. 'Oil ). In fact, 
stars of 0(100) Mq have very short lifetimes and either 
explode as pair instability Sup er Novae (PISN e) or di- 
rectly collapse to black holes (jHeger et al. '03[ ). They 
also contribute to the metallicity enrichment of the early 
universe, shaping the transition to the second genera- 
tion of stars. Moreover, the very high temperature of 
these objects makes them efficient engines for the pro- 
duction of Lyman- Werner ultraviolet photons, thus ini- 
tiating the cosmic rcionization process. Unfortunately, 
all of the traces PopIII stars leave behind are challeng- 
ing to observe and so far no unambiguous detection has 
been recognized, although the detection of anisotropics 
in the infrared background consistent with the exi stence 
of PopIII has been claimed ([Kashlinskv et al. '05f ). 

Despite the difficulty of detecting neutrinos, they may 
be potentially interesting messengers of the high-redshift 
universe since, differently from gamma rays, neutrinos 
travel unimpeded over cosmological distances. Also, 
the peculiar initial mass function (IMF) of PopIII stars 
should increase their Super Nova (SN) rate with respect 
to later stellar populations. The peculiarity of PopIII 
stars may well produce high energy neutrinos above the 
expectations for PopII stars, thus making them an in- 
triguing target for neutrino telescopes. A fairly general 
expectation is that PopIII may emit large amounts of 
neutrinos , both thermally produced at the time of the 
collapse (llocco et al. 'OSf i and non-thermally during a 
Gamma -Ray Burst (GRB) p hase associated with the ex- 
plosion ([Schneider et al. '02l ). GRBs are one of the can- 
didate hadronic accelerators, and any proton accelerator 



is a potential emitter of high energy neutrinos, produced 
via hadronic (pp) or photo-hadronic (^7) reactions in the 
surrounding medium. Although the details of this whole 
picture are far from being established and despite the 
large uncertainties in the models, the idea of using neu- 
trinos to probe PopIII stars seems promising, and surely 
deserves further study. 

In this paper we perform a new estimate of the high en- 
ergy neutrinos expected from GRBs associated with the 
first generation of stars in light of new models and con- 
straints on the epoch of rcionization (EoR) and a more 
detailed evaluation of the neutrino emission yields. The 
goal of this paper is two-fold: {i) to compare the diffuse 
high energy neutrino background from PopIII with the 
one from "ordinary" (PopII) stars that contribute with 
PopIII to the rcionization at 7 < z < 12; {ii) to discuss 
the chances of detection on the light of the performances 
expected for current or future neutrino telescopes as a 
function of the astrophysical input. Indeed both con- 
ditions need to be fulfilled to establish if there is any 
realistic chance to discriminate the high energy neutrino 
emission from PopIII. We anticipate that we find that 
high energy neutrinos from Population HI stars will not 
be observable at current or near future neutrino tele- 
scopes, falling below both IceCube sensitivity and atmo- 
spheric neutrino background under the most optimistic 
assumptions for the GRB rate. The disagreement with 
previous literature is mostly due to unrealistically high 
PopIII GRB rates previously considered, following e.g. 
from values of the ionization efficiency nowadays con- 
sidered too extreme based on self-consistent rcionization 
scenarios. 

The plan of this paper is the following: after introduc- 
ing the basic formalism in Sec. [21 we shall preliminarily 
discuss what a semi-empirical estimate of the GRB rate 
at high redshift can tell us (Sec. [3|). Then, we devote Sec. 
[3|to describe the models used to perform theoretical esti- 
mates of the neutrino fluxes. We shall pay particular at- 
tention to compare estimates derived consistently within 
the same cosmological assumptions. In Sec. [5|we present 
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our results and compare with the chances of detection, 
and finally in Sec. [6] we conclude. In Appendix [A] we 
report some detail of the GRB model used to compute 
the neutrino yields. 

2. BASIC FORMALISM 

In this Section we establish the formalism we will use 
throughout the paper to estimate the diffuse flux of neu- 
trinos emitted by GRBs. The integrated signal observed 
today at energy Ei, is 

£;2$^ [GeVcm-2g-ig^-i] dF. 



'dE, 



dzdn dE^dA' ^'^'^ 

where dNc/dzd^l is the differential rate of GRBs which 
beam towards us per unit solid angle, and dNl^°/dEiydA 
is the average flux emitted by a single source at energy 
E^{1 + z). The function dNQ / dzd^l can be written as 

b G{z) dV 



dNG 1 pg{z) dV 

:v-) — 



(2) 



dzdl^'"' An{l+z)dz 47r (1 + z) dz 

where Pg{z) is the GRB rate in a comoving volume and 
dV is the comoving volume element, such that 



d7 = ^"if(I)^(^)' 



(3) 



We have introduced the comoving distance r[z) defined 
as 

c 



r{z) 



■ dw . 



(4) 



/o H{w) 

The Hubble function H[z) writes in terms of the frac- 
tions of the critical energy density in matter and 
cosmological constant Da as 



H{z) = Ho^Um^I + zf + + z)2 + r!A , (5) 
Hq being the inverse Hubble distance today H^^ ~ 
3 /i^^Gpc and fifc = 1 — VLm — ^k- We have also denoted 
by 6 = (1— cos6'jot) the averaged beaming factor of a jet of 
opening angle so that G{z) is the beaming-corrected 
overall GRB rate. The quantity h also fixes the ratio be- 
tween the physical energy Ejot and the one released by 
the GRB if that was isotropic, 6/47r = Ejct/iJiso- If we 
introduce the function of the energy at the source E'^j 

Td7V^° 



dE' 



(6) 



the average energy flux, El{dNl^° / dE^dA) (units GeV 
cm~2), observed at a given energy Ei, and emitted by a 
source at redshift z can be expressed as 



E 



dN'J° 
" dE,AA 



1 J„[F,(l + z)] 



(7) 



l + z AT:r^{z) 

In the previous expression, E^{1 + z) is the emission en- 
ergy, Jy[Ei,{l -\- z)]/{l 4- z) takes into account the red- 
shifted energy spectrum, and possible boosting factors 
are already included in J^. 
Putting the two pieces together, we finally get 

, _cb f ME{l + z)]G{z) 

^•^^^ -4^ J {i + z)m{z) ■ 

It remains to estimate G{z) and J^{E), a problem which 
we shall address in the following. 



3. AN EMPIRICAL ESTIMATE OF GRB RATE 

Let us start by using an empirical approach: the 
GRB rate is derived from o bservations and ext rapo- 
lated at high redshifts, as in (fYonetoku et al. '031 ) and 



(jFenimore fc Ramirez- Ruiz '05|). This approac h assumes 
the E'p-luminosity relation, ( Amati et al. '02l ). to infer 
the redshift of a GRB whose host is at unknown redshift. 
Note that in order to extrapolate the i?p -luminosity re- 
lation to high redshifts, one implicitly assumes that al- 
though GRBs at high redshift might belong to different 
stellar populations with different star formation rates 
(SFRs) and IMFs, they constitute "standard candles" 
carrying no memory of the population of the progenitors 
they belong to. 

Aware of this caveat, the most natural hint of a transi- 
tion PopIII-PopII would then be a break of the power-law 
dependence of the GRB formation history on the red- 
shift; indeed, one should expect in general that the con- 
ditions for the GRB appearance will occur at a different 
rate in stellar populations with different characteristics. 
This is manifestly not the case of t he GRB rate derived 
empi rically from observations in ( Yonetoku et al. '03l ) 



and (jFenimore fc Ramirez- Ruiz '05 ): no discontinuity 
appears at 6< z <12, where a PopIII-PopII transition 
would suggest to see it. In particular, these results in- 
dicate that the GRB formation rate always increases to- 
ward z ^ 12, in a way parameterized by: 

/ ^ _ /(I + ^)": z<l 

Pcmp[Z) — Pcmp[y}) X "j 2"-/3(l + z)'^ 1 < Z < 12 ^ ' 

where a = 6.0 ± 1.4, /3 = 0.4 ± 0.2 and 

Pomp(O) - 2.3 X lO-i-'s-iMpc-^ , (10) 

all the quantities being fixed according to 
(jMurakami et al. '05f ). where for the normalization 
we use their intermediate estimate for the luminosity 
correction {k = 0.5). Note that if we assume a simple 
top- hat function for J^,, from Eq. ([5]) and the fit for 
Pcmpiz) it is easy to check that neutrinos emitted at 
high redshift have at most a subleading contribution 
to the total flux. Q ualitatively, this is consi stent with 
what found e.g. in (jMurase fc Nagataki '05[ ): most of 
the GRB neutrino signal is not sensitive to the GRB 
rate at high redshift. 

This general result might be interpreted as (weak) 
empirical evidence for the dominance of PopII progen- 
itors in the overall GRB population at any redshift. 
However, this conclusion only holds under two assump- 
tions: (i) the validity of the iJp-luminosity relation up 
to high redshifts; since this is an empirical relation and 
the engine of GRB is far from being understood, there 
is no guarantee that this is the case; (ii) the absence 
of a significant fraction of choked GRBs among PopIII 
stars. It has been argued that due to the different struc- 
ture and properties of PopIII stars the jet in a coUap- 
sar may be una ble in most cases to p unch through the 
stellar envelope (jMacFadven et al. 'Oil ), thus choking the 
burst in gamma rays, although still producing large neu- 
trino yields. This w as also the working hypothesis of 
(jSchneider et al. '0^ . These arguments suggest that a 
more realistic estimate of the neutrino flux from PopIII 
GRB may require some degree of modeling, which we are 
going to address in the next section. However, we can 
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anticipate that the empirical GRB rate obtained in this 
section is roughly in agreement with the one we obtain 
from the CF05 model we shall introduce in the following 
section. 



4. THEORETICAL ESTIMATES OF THE NEUTRINO FLUX 

Let us now estimate the GRB rate at high redshifts 
from current theoretical models on the stellar popula- 
tions at high redshift, under the assumption that the 
GRB rate tracks the SFR. Besides the SFR for the 
different stellar populations, it is clearly mandatory to 
know their IMFs in order to evaluate the fraction of 
stars that are likely to end their lives as coUapsars, 
believed to give rise to a GRB. In this section we 
present the neutrino fluxes expected from the two dif- 
ferent populations of stars under different EoR models 
available in literature. The three models of reionization 
we consider are su mmarized in Table [H the "fiducial" 
model presented in (jChoudhurv fc Ferrara '05l ). figure 1 
at page 586 (hereafter CFOSq); the one in figure 4, page 
590 of the same paper (here after CF05{,); the model 
in (jChoudhurv fc Ferrara '06h (hereafter CF06), which 
takes into account the new cosmological data on reion- 
izatio n from the third-y ear data release of the WMAP 
team (jSpergel et al. '06[ ). These models are used in the 
following to derive the star formation rates needed for our 
estimates. All the CF-models share the same physical 
assumptions: three radiation sources are taken into ac- 
count during EoR: Quasi Stellar Objects (QSOs), PopII, 
and PopIII stars which contribute with different char- 
acteristics to the radiative and chemical enrichment of 
the IGM, and to the consequent feedback. Among the 
many parameters of the models, a key role is played by 
the fraction of escaping photons per halo (directly re- 
lated to the ionizing efficiency), which must be deduced 
from the knowledge of the source and halo characteris- 
tics of PopII and PopIII stars; the luminosity function of 
QSOs is calculated according to standard assumptions. 
Both the CF05 models assume that all PopIII stars have 
a mass of 300Mq while PopII follow a Salpeter IMF, 
S{m). It is also worth mentioning that the transition 
between PopIII and PopII stars, regulated by the chemi- 
cal feedback, is considered to be instantaneous. The only 
difference between CF05a ad CF05f, is the ionizing effi- 
ciency of PopIII, which in the latter model is assumed 
to be 2/7 of the former one. This drastically changes 
the normalization and shape of the PopIII SFR. With 
respect to the CF05 models, in CF06 the stellar forma- 
tion in minihaloes is suppressed, and a self-consistent 
calculation is made of the chemical enrichment process, 
which takes longer and results in a non-instantaneous 
transition between PopII and PopIII. For our purposes, 
however, the crucial difference is that CF06 assumes a 
Salpeter IMF for the PopIII, too. The authors imple- 
ment this choice in order to match the NICMOS high 
redshift source counts. We address the reader to the 
original papers for more details. Also note that none of 
the models provides the low-redshift (z < 3) extrapola- 
tion of the SFR. Therefore, for the neutrino fluxes from 
PopII shown as comparison in the following, we have 
used at low-redshift the GRB rate given in Eq. (12c) of 
(H urase fc Nagataki '05[ ) normalized in order to match 
the observe GRB rate z=0. 



Model 


Refcrenee 


Notes 


CF05a 


(Choudhurv H'. Ferrara '05) 


high ion. efficiency 


CFOSfc 


(Choudhurv & Ferrara '05) 


low ion. efficiency 


CF06 


(Choudhurv & Ferrara '06) 


WMAP3, 5(M) 



TABLE 1 

The models of reionization we consider for calculating 
the neutrino fluxes. see text for details. 



4.1. Rate of GRBs from PopII 

If we denote by i?n the SFR of PopII stars (taken from 
the models reported in Table[T]), an estimate for the num- 
ber density rate of GRBs from PopII stars, Gn, can be 
obtained as 

Gu{z) = anlu^ (11) 
Mil 

where AIu is the average mass of a PopII star, an is the 
fraction of PopII which are likely to form core collapse 
SNe and 711 is the fraction of core collapse SNe which are 
likely to form a GRB. The average mass of a PopII star 
is obtained via an integral over their IMF, Tu, as 



A/ti 



/ Xii(A/)A/dM 
/Xii(A/)dA/ 



(12) 



assuming (consistently with all the EoR models con- 
sidered) that Tii(Af) is a Salpeter mass function 5'(A/) 
which writes (jKroupa '00[ ) 



S{m) oc m 



0.3, 0.01 <TO< 0.08 

1.3, 0.08<TO<0.50 (13) 

2.3, 0.50 < TO, 



where m is the mass in solar units A/q. This is also used 
to estimate an, assuming that all PopII with to > 10 will 
end their lives as SN, 



CTii 



J^^'^ S{m)dm 



(14) 



We allow for a z— dependence in 711 to take into ac- 
count the GRB-metallicity anti-correlation and the fact 
that the m etallicity evolves with the a ge of the Univer se. 
Following (|Yuksel fc Kistler '06D and (|Yoon et al."^ . 



7ii(^) = 7ii(0)(l + z)' 



(1 



\1.4 



1250 



(15) 



where we fixed 711 (0) according to (|Yoon et al. '06f ) and 
7ii(z) has been considered constant for z > 4.5, thus 
considering no metallicity evolution above this redshift 
(however, this choice affects very little the final results, 
most of the contributions to the GRB neutrino back- 
ground coming from bursts that occur ai z < 7). 

4.2. Rate of GRB from PopIII 

Formally, the number density rate of GRB from PopIII 
stars, Giii can be written identically to Eq. (jlip . namely 

Riuiz) 



Giii(z) = criii7in ■ 



Mm 



(16) 



where the various coefficients and functions involved 
now assume different values. The Riii{z) are taken 
from the models reported in Table [TJ Concerning the 
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Model 


Mass (Mq) 


Eiso (erg) 


Ejct (erg) 


r'is (cm) 


'1 


Hiddcn-A 


200 


lO''-* 


1052 


1013.0 


1 


HIddcn-B 


200 




lO^i 


10i3.u 


1 


Hiddcn-C 


60 




10^1 


1012.5 


1 


Hiddcn-D 


60 


1053 


IqSO 


1012.5 


1 


Hiddcn-E 


200 


10'>4 


10^2 


lOii.o 


10 


Hiddcn-F 


60 


1053 


10^1 


1011.5 


10 


Prompt 


35-125 


1053 


1.24 X 10^1 


1013-15.5 


1 



TABLE 2 

The parameters varied in the models of GRB considered. 



IMF of PopIII stars, Tni, there arc currently two hy- 
potheses: (i) a bimodal IMF with a first peak at a 
few solar masses and the second one at ^ 100 Mq 
(|Nakamura fc Umemura 'Oil ): (ii) a peaked IMF around 
the range of few hundred solar masses has been later pre- 
sente d in (lAbel et al. 'Olfl . Recently it has been pointed 
out, (jO'Shea et al. '05D . that primordial stars forming in 
a halo neighboring the formation star of a PopIII have 
still zero metallicity but smaller accretion rates onto the 
central dense region, resulting in a final stellar mass of 
10-100Mq. This "intermediate" generation would have 
equally contributed to the reionization and we do im- 
plement this hypothesis by studying the contribution of 
M ~ QO Mq in our neutrino models. Whatever the de- 
tails are, the hypothesis (ii) is nowadays widely accepted 
in the community, or in any case it is the Xni most fre- 
quently implemented in EoR models. We shall consider 
the following cases: 



(a) 



Jiiioc(5(Miii-3OOAf0), 

i.e. a toy model of a "monochromatic" high mass 
mode for the IMF. 



Ziii(X(5(A/iii-60Mq), 

i.e. a toy model of a monochromatic IMF, with the 
average mass of metal-free stars forming from the 
collapse of halos triggered by the explosion of the 
very first stars. 



Xiii (X S{M) , 

as assumed e.g. in (jChoudhurv fc Ferrara '06h . 

Note that in both case (a) and case (6) one has ctih ~ 1, 
while in case (c) one obtains an expression analogous to 

Eq. mi). 

Concerning 7111, there arc no firm estimates. For the 
sake of clarity and to simplify the c omparison with previ- 
ous results of (jSchneider et al. '02| ). we shall assume the 
same value, 

7111 ^ 1 . (17) 

This corresponds to an upper limit to the PopIII GRBs; 
a smaller GRB efficiency would result in a down-scaling 
of the results presented in the following. 



4.3. Neutrino yields from GRBs 

High energy neutrinos from GRBs are predicted in var- 
ious scenarios. The one most frequently discussed is 
neutrino emission in the internal shocks that produce 
prompt emission ([Waxman fc Bahcall '97f ). In this sce- 
nario, accelerated protons interact with gamma-rays via 
photomeson production, and produce pions and kaons, 
which decay into neutrinos. The observed gamma- 
ray spectra are usually represented by a broken power- 
law. As for proton spectra, a first order Fermi- type 
spectrum is frequently assumed. This possibility has 
been studied by many a uthors. Murase fc Nag ataki 
(jMurase fc Nagataki 'Ol iMurase fc Nagataki '06[ ) did 
such calcultions using Geant4 with experimental data. 
They also took into account various cooling processes. 
One parameter set o f their results (which is used in 
(|Achterberg et al. '07D ) will be shown for comparison in 
this paper, and named Prompt in the following. 

Anothe r model was first suggeste d by Meszaros and 
W axman (iMeszaros fc Waxman 'Olh and later extended 
in (jRazzaque et al. '03f ). The GRB progenitor is usually 
taken to be a massive star with a He core and H envelope. 
The central engine of GRBs is still a major open prob- 
lem in astrophysics; the leading model for (long duration) 
GRBs is a huge stellar collapse leading to formation of a 
highly rotating black hole with an accretion disk (coUap- 
sar model). The jet which produces succeeded GRBs or 
choked GRBs propagate in stars. This jet would inter- 
act with the stellar envelope before it produces prompt 
emission. In addition, the observed submillisecond vari- 
ability allows us to expect that intenal shocks can occur 
at sufficiently small radii, where the Thomson optical 
depth exceeds unity. Protons can be accelerated in such 
internal shocks, and they can interact with photons from 
electrons that are accelerated in internal shocks and/or 
termination shocks^ . Photons from such inner radii can- 
not escape due to large optical thickness, so that they are 
hidden sources as far as gamma-rays are concerned. Only 
neutrinos would be useful as the probe of physical pro- 
cesses in this region. Since here we want to estimate the 
high energy neutrinos from PopIII GRBs, consistently 
with the empirical suggestion from Sec. [3] we shall only 
evaluate neutrino emission from such inner radii. In- 
deed, this emission is present also for chocked GRBs, 
and we stress that such "hidden" models are the ones 
that best match our choice 7111 ~1. In fact, models such 
as the Prompt would predict a lower ratio of GRBs, thus 
creating friction with our choice. Throughout the pa- 
per, we take into account neutrino oscillations in vacuum 
by adopting mixing angles 612 = 0.59 , 623 = 7r/4 , and 
013 = 0. The yields shown in the rest of this paper refer 
to the muon neutrino flavor, which is the best channel for 
detection of a low-statistics signal at neutrino telescopes 
given the long tracks of muons produced in charged cur- 
rent interactions. For diffuse signals competing with the 
atmospheric background, for which the separation must 
be done on the basis of the different energy spectrum, 
the electr on neutrino showers m ay be an interesting ob- 
servable (jBeacom fc Gandia '04[ ). However, showers are 
only detected if contained in the instrumented volume, 
which lowers the effective collecting area of the telescope 

^ More details about the terminology and GRB modeling are 
given in Appendix lAl 
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with respect to the muon track case. Wc shall see that 
our event rate predictions for muon tracks fall below the 
sensitivity of a km^ instrument, which makes hopeless 
exploring alternative channels. 




Fig. 1. — Ji, for the different models used, see Tab. [2] 



Since our knowledge of the GRB engine and 
PopIII properties is still limited, it is unclear 
how efficient are PopllI stars as GRB progeni- 
tors. But rotation and chemical mixing seem to 
be a key ingredient towards successful co Upsar mod- 
els (jMaeder fc Mevnet 'OOl: iHeger et al. 'OOO . and recent 
simulations in ( Yoon et al. '06D suggest that GRBs can 
be born from PopIII stars, an hypothesis we will assume 
throughout. We expect that m ore massive star collapses 
lead to larger released energies (|Heger et al. '03f ). Obser- 
vationally, to ease the energy requirements, it is usually 
considered that a globally asymmetric, relativistic jet is 
launched from GRB progenitors. In fac t, this picture i s 
suggested by some observ ations, e.g. (iFrail et al. 'Olh . 
The Ghirlanda relation (jGhirlanda et al. '04l ) implies 
that bursts can have various jet energies, and some 
GRBs may have Ejot ^ 10^^ erg, which is larger than 
the frequently used value E-^ct ~ 1-24 x 10^^ erg. The 
isotropic energy is highly uncertain, too, and it may 



be as high as Ei^ 



10 erg; we adopt this value 



also to allow direct c omparison with previous results in 
(jSchneider et al. '02f ). 



Because of these high uncertainties, we explore a wide 
range of parameter values for E-^^t and -Eiso, as reported in 
Table [21 The other parameters varied are 77, a factor de- 
termining the maximum acceleration energy of protons, 
and Tig. Further details on the model and its parameters 
are in Appendix [A] 

The radial density profiles 5{r) for the PopIII pro- 
genitor mass-models, M = 200Mp, and M = e OMp,, 
have been calculated in (jHeger fc Wooslev '02l ) and 
() Rockefeller et al. '06( ). respectively. We would like to 
comment on the apparent contradiction between the 
Iiu oc S{Miii — 300 Mq) motivated in the previous sec- 
tion and the M = 200 Mq model used for the calcula- 
tion for the neutrino flux. Physically, stars with masses 
M = 140 Mq <M< M = 260 Mq arc known to directly 
explode as PISNe without the formation of a Black Hole, 
thus being inefficient GRB progenitors. While important 
for the GRB rate, however, we expect that the differ- 
ences in the density profile for the two star masses are 
much less pronounced in the neutrino yields, and cer- 
tainly within the uncertainties and assumptions made in 
such an uneven ground. 

In Fig. [1] we show the neutrino emission spectra at 
the rcdshift of the source Jiyjii{E) from the PopIII GRB 
models considered in Tabic [2l also compared with the 
spectrum Ji/ji(£') for the PopII model Prompt. 

In the scenario where internal shocks and termination 
shocks occur at sufficiently small radii inside the pro- 
genitor star, almost accelerated protons are depleted by 
hadronic and photomeson reactions. However, produced 
mesons and muons suffer from cooling processes. The 
strong magnetic field and copious photon field make pi- 
ous and muons cooler before they decay. Hence, the neu- 
trino flux will be suppressed due to cooling processes. In 
Fig. [21 the relative fiux level basically reflects the differ- 
ences among isotropic energies of models. However, the 
strong field makes mesons and muons give up their en- 
ergies, so that such suppression becomes important. For 
example, the difference on the flux level between Model 
E and Model F is reduced. Furthermore, the strong field 
can force the higher break energy, which is determined 
by equating the pion lifetime and its cooling time, to 
be shifted to lower energies. This effect also leads to a 
double peak structure as shown in Fig.l for Model A- 
D. The lower peak corresponds to the contribution from 
the termination shock. Here, in fact, nonthermal protons 
accelerated in the internal shocks interact with thermal 
photons produced by electrons accelerated in the termi- 
nation shock. On the other hand, the higher peak id due 
to protons from the internal shocks themselves. If we 
assume similar equipartion parameters, the termination 
shock can lead to the stronger field than that of inter- 
nal shocks. Therefore, neutrinos from this region have 
lower energies due to pion and muon cooling and this 
explains the double peak structure. However, the promi- 
nence of this feature is probably due to the simplicity of 
the model considered: in realistic situations, the emission 
regions would not be approximated by the simplest two 
zone model (internal shock region and termination shock 
region) , and the gradient of the field strength should be 
taken into account. Then, the double peak structure 
would be more smoothed, although a detailed treatment 
of this issue goes beyond the scope of this paper. 
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Fig. 2. — The E^^^^m for different J^^m models, assuming 
CF05a as EoR model. 

5. RESULTS AND DISCUSSION 

In this section we present our results for the neutrino 
fluxes expected at the Earth. In Fig. [2] we show the 
contribution of PopIII stars calculated as described in 
the previous sections, assuming CF05a as a fiducial EoR 
model. In general, these fluxes present a huge variability 
due to the large uncertainties on the production parame- 
ters, yet they are all much lower than the most optimistic 
results found in (jSchneider et al. '0^ . 

In Fig. [3] we compare the high energy neutrino flux 
at Earth expected from PopII stars with the one from 
PopIII for different but mutually consistent EoR mod- 
els, as denoted by the same line-style. When accounting 
for the assumed baryon loading factor of ^acc = 10 (see 
Appendix for the role of this parameter) the estimated 
contribution from PopII stars agrees with typical results 
found in the literature, and the dependence on the EoR 
model used is marginal, reflecting the fact that at low 
redshifts (where most of the signal comes from) all mod- 
els must agree with the available data ^. On the other 
hand, the situation for PopIII is different. The solid and 
dashed lines assume the model A neutrino yields, which 
as shown in Fig. [2] maximizes the neutrino production. 
The dotted line instead follows from the CF06 model, 

^ The slight discrepancy with previous results in the PopII neu- 
trino flux at Earth is due to the higher GRB rate we use for 2 > 3, 
arising from the reionization models used here. 



where we have implemented a Prompt spectrum for J^ju 
coherently with the assumption of a Salpeter IMF. We 
have treated this case as the PopII one, thus obtaining 
the fraction of coUapsars 7ni — 10~^, consistently with 
no mctallicity evolution beyond z = 4.5. This illustrates 
the high energy neutrino flux for the case of an hypothet- 
ical low-mass PopIII generation; whereas the CF06 SFR 
coupled with the A model for J^ju gives results almost 
indistinguishable from CF05f, and we have therefore not 
plotted it. The solid and dashed lines show that the SFR 
has only a minor effect on the neutrino flux at Earth; on 
the other hand, their comparison with the dotted line 
shows that the IMF plays a fundamental role in shap- 
ing the ffux at Earth, affecting the fraction of stars that 
will give rise to a collapsar and the magnitude of their 
explosion energy. 

A fairly robust outcome of our study is however that 
neutrinos from PopIII GRBs will not be detectable with 
current or near future experiments. In fact, the PopIII 
flux shown in Fig. [3] falls below both the cu rrent 
AMANDA-II bound (4-ycars data, see (jHalzen '06l )) and 
future prospects for five years of IceCube exposure. Even 
worse, this contribution is buried beneath the atmo- 
spheric neutrinos, whose average spectrum and uncer- 
tainty is plotted in F ig [2l according to the compilation 
recently reported in (|Evoli et al '07f ). Also, it is worth 
stressing that while the PopII fiux we have plotted in 
Fig [3] is a relatively robust expectation because the back- 
ground can be evaluated from the observed GRB rate, 
the PopIII one is an upper limit, as it assumes that 
all the stellar events will give rise to a collapsar with 
high energy neutrino emission. These results greatly dis- 
agr ee with the much mo re optimistic estimates obtained 
in ([Schneider et al. '02t ). thus motivating some explana- 
tion. We believe that our discrepancy with previous re- 
sults is to impute to the different no rmalization of PopII I 
SFR. In particular, from Fig. 1 in ([Schneider et al. '02l ) 
one can infer that the neutrino spectra results plotted in 
their Fig. 2 require a rate of ~ 10^ GRBs per year versus 
e.g. a rate of ~ 10^ GRBs per year for the C FOSg mod el 
adopted here (more details can be found in ([locco '07l )). 
Actually, such a large discrepanc y only arises for th e 
most extreme model considered in ([Schneider et al. '02l ). 
which adopts a value for the ionizing efficiency which is 
far too large, based on the results of self-consistent mod- 
eling of the reionization epoch nowadays available. Pre- 
dictions for more realistic models are in better agreement 
with our results. It is also interesting to note that the 
theoretical models used here are consistent, within a fac- 
tor of a few, with the estimate for the overall GRB rate 
obtained empirically in Sec. |3l which makes us confi- 
dent on the more realistic estimates provided here with 
respect to previous literature. 

6. CONCLUSION 

We have performed a study of the high energy neu- 
trino diffuse background which is to be expected from 
Population III stars, under the assumption that they 
will end their lives as GRB. We have compared it with 
the analogue expected from Population II stars, using 
mutually consistent SFRs obtained from EoR models 
available in literature. Our estimate of the PopII GRB 
rate has been performed under widely accepted assump- 
tions and leads to estimates of the high energy ncu- 
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Fig. 3. — The muon neutrino flux for the different Reion- 

ization models, using the A model for J^jn. Upper-right set of 
curves refer to the PopII flux, lower-left ones to the PopIII contri- 
bution under optimistic assumptions (see text). 

trino fluxes which are in agreement with previou s re- 
sults (jMurase fc Nagataki '051 : [Razzaaue et al. '03f ). On 
the other hand, we have presented a maximal model for 
PopIII stars, assuming that all of them will end their 
lives with a GRB, either choked or not. Even under such 
optimistic assumptions, a detection of the diffuse high 
energy neutrino background expected from Population 
III stars appears out of reach. Even worse, the contri- 
bution from PopII GRBs would contaminate that from 
PopIII GRBs. In addition, we cannot expect neutrino 
signals correlated with gamma-rays from GRBs, because 



it is thought to be very rare to sec PopIII GRBs by the 
current satellite such as Swift. Hence, the PopIII neu- 
trino signals are expected to be hidden by atmospheric 
neutrino background. With reasonable values of the non- 
thermal baryon loading factor, an extreme contribution 
from PopIII GRBs — including possible choked bursts — 
falls in fact underneath IccCube five years sensitivity 
and is overwhelmed, at low energies, by the atmospheric 
neutrino background. This implies that although PopIII 
stars may contribute to the high-energy diffuse neutrino 
background and their spectrum would be indeed sensi- 
tive to their IMF and SFR, neutrinos cannot be used as 
a diagnostic tool to check properties of either the pop- 
ulation of stars during the epoch of Reionization or the 
GRB internal shock properties, thus confirming the elu- 
sive nature of the earliest generation of stars. 



ACKNOWLEDGMENTS 

F.I. thanks A. Heger and C.L. Fryer for kindly pro- 
viding the stellar models, and A. Ferrara for kind clari- 
fications about EoR models. P.S. acknowledges support 
by the US Department of Energy and by NASA grant 
NAG5-10842. This work is in part supported by a Grant- 
in-Aid for the 21st Century COE "Center for Diversity 
and Universality in Physics" from the Ministry of Educa- 
tion, Culture, Sports, Science and Technology of Japan. 
S.N. is partially supported by Grants-in-Aid for Scien- 
tific Research from the Ministry of Education, Culture, 
Sports, Science and Technology of Japan through No. 
19104006 and 19740139. K.M. is partially supported by 
Japan Society for the Promotion of Science. 



APPENDIX 

MODELING NEUTRINO EMISSION FROM JETS INSIDE THE GRB PROGENITOR STAR 

Here we summarize the basic mechanisms by which neutrinos are produced in jets inside the GRB progenitor star. 

We assume that a (long) GRB is a huge stellar collapse leading to the formation of a highly rotating black hole with 
an accretion disk, during which a globally asymmetric, relativistic jet is launched from the GRB progenitor. The GRB 
observations imply that the Lorentz factor of jets is very large (Fjet ~ 100) in the gamma-ray emitting region. This 
is the final value which is achieved outside the star, and it may be representative of the intrinsic injection Lorentz 
factor. Of course, we do not know the intrinsic value and it may also be much smaller than the final value. Wc shall 
briefly comment on the latter case later. When the jet propagates inside the star, it will produce a bow shock ahead 
of it. This jet is capped by a termination shock and a reverse shock. The jet termination radius (where the jet is 
decelerated) is indicated by rh. If the jet is highly variable and the variability time scale St is small enough, internal 
shocks can also occur inside the pre-decelerated jet. The internal shock radius is written as r^s ~ 2T?^^.cSt. In the 
usual internal shock scenario of GRBs, the observed gamma-rays are attributed to internal shocks which occur in the 
optically thin region outside the progenitor star. Lower Fjot and/or St lead to smaller radii below the stellar surface, 
which will be optically thick. Here, we are interested in such opaque, sub-surface internal shocks, and we assume 
fis < f*, where is the radius of the progenitor's stellar surface. 

These internal shocks are expected to be coUisionless, so that both electrons and protons may be accelerated. If the 
magnetic field is strong enough, electrons can be accelerated up to high energies and radiate photons. The photon 
density in the pre-decelerated jet is given by 

where Al is the width of subshcUs, for which wc use Al w r/Fjot because we have assumed that the jet acceleration 
has already ceased. The magnetic energy density in the field B is expressed as 

C^i'' = f^=ei°*f/f . (A2) 

Due to the existence of the strong magnetic field, electrons radiate synchrotron photons. These photons will be 
thermalized; if the jet is optically thick to Thomson scattering. We assume = 0.1 and expect that the jet is 
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opaque. Therefore, we can approximate the spectral distribution of the radiation of energy density Ui^"^ by a black- 
body spectrum, and associate to it a temperature Tjct- The above approximate treatment is accurate enough for our 
purpose. 

Protons get accelerated in the shocks as well. It is widely believed that cosmic rays can be accelerated by the first- 
order Fermi acceleration mechanism. We assume that this mechanism can work efficiently and adopt the spectral index 
~ 2. The non-thermal proton energy density is expressed as {7^,°* = ^.accUi^^*", where ^acc is the nonthermal baryon 
loading factor and for efhcient proton acceleration we can express as ^acc ~ l/ce, Ce being the fraction of internal 
energy carried by the electrons. Because plausible values of this parameter is not known yet, we adopt ^acc = 10, 
which correspond s to the assumption that the energy of protons pe r logarithmic energy bin is comparable to the GRB 
radiation energy (jWaxman fc Bahcall '97l ; [Murase fc Nagataki '05[ ). Too large value of face requires too small values 
of ee, which is usually unexpected in GRBs (but there is no proof because we do not know the total explosion energy). 
The non-thermal proton spectrum is given by 

" p = P c-^ (A3) 

dsp ln(e™a'^/e™'") ^ ^ ' 

The minimum energy e™'" is set to 10 GeV. The choice of this value is not so sensitive to our final results. On the 
other hand; the maximum energy is important for the purpose of knowing neutrino spectra at the highest energies. 
The maximum energy is determined by the condition 



eBc 



(A4) 



where at the l.h.s. it appears the Larmor radius of the proton times the pre-factor r]{= 0(1 — 10)) which depends on 
the details of acceleration mechanism; taee is essentially the acceleration time in the Bohm limit, and the sum at the 
r.h.s. extends over all the energy loss channels timescales U. 

Accelerated protons can interact with protons themselves or photons, which lead to meson production such as 
pions and kaons. In this pape r, we consider neutrinos produced th rough photomcson production only. The method 
of calculation is described in (|Murase fc Nagataki '05l : iMurase 'Q7( ). Here, we briefly sketch this process by simple 
analytic considerations. The photomeson production is a threshold process, with a threshold energy of about 145 
MeV in the rest frame of the incident proton. The dominant inelastic channel is P7 A+ with the cross section 

10 cm and the inelasticity Kp w 0.2, which is called A-resonance. For the cases we consider, the 
optical depth at the A-resonance is very large, so that the photomeson production efficiency can also be high. This 
means that almost all the protons that have sufficiently high energies (above the threshold energy) will be depleted 
due to photomeson production. 

The threshold for pp inelastic interactions is lower, so these processes will also occur. Although this may be an 
important neutrino source in the TeV energies (in particular for low Lorentz factors Fjot), in this paper we focus 
on high energy neutrinos produced by sufficiently high energy protons above the threshold energy for photomeson 
production. We take into account the pp process only for estimating the proton energy loss time scale, treating this 
process analytically for simplicity. 

Next, we consider the interaction between the jet and the progenitor star. As the jet advances through the star, it 
drives a bow shock ahead of it. The jet is capped by a forward shock, and a reverse shock moves back into the jet, 
where the relativistic jet is decelerated. The shocked jet plasma and shocked stellar plasma would advance together 
with a jet head Lorentz factor Fh Tjet- By equating the pressures behind the forward shock and reverse shocks, one 
finds the following estimate for the Lorentz factor of the jet head, 

\ 

(A5) 




where mp is the proton mass and S is the mass density of the environment. 
The relative Lorentz factor between the shocked jet plasma and un-shocked jet plasma is 

1 / Tjct , Fh 



Electrons will be accelerated in these shocks. However, the electrons would give up all their energy on a very 
short time scale, by synchrotron and inverse-compton (IC) cooling, converting a large fraction of the shocked plasma 
internal energy into radiation. These radiated photons in the shocked jet plasma will be thermalized due to large 
optical thickness. Hence, target photon density will be approximately a black-body radiation, with an overall energy 
density 

U!; ~ (4F' + 3)(F' - l)nf mp , (A7) 

with an associated temperature Th. 

The reverse shock is likely to become radiation dominated, so that the IC cooling by the electrons becomes important 
and affects the dissipation of jet kinetic energy. If the reverse shock is indeed radiation dominated, the shock thickness 
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would also be of order the mean-free path of thermal photons propagating into the jet. The copious photon field also 
affects the neutrino spectrum. Charged mesons produced via P7 (and pp) interactions will suffer IC and synchrotron 
losses, and neutrino spectra will be suppressed. The magnetic energy density in the shocked jet frame is expressed as 

Uh=eBU^- (A8) 

For = 0(0.1 — 1), the magnetic field strength is very large, up to ~ 10^~^ G. Hence, charged particles would 
suffer from synchroron loss, and one should take into account synchrotron losses of pions in order to calculate neutrino 
spectra, as we do in our computation. 

If non-thermal protons that are accelerated in internal shocks are not completely depleted, they can enter the 
shocked jet. Such protons are expected to interact with photons in the shocked jet region, and to produce neutrinos. 
In principle, protons can be accelerated at the termination shocks, too. However, cooling processes such as IC loss 
significantly will reduce the maximum proton energy there. Hence, we treat the contribution from protons that are 
accelerated in internal shocks only. The proton spectrum in the frame of the shocked jet plasma is given by 

where /j,''' is a fraction of depleted protons in the internal shocks, which can be estimated as the ratio of the dynamical 
timescale involved with the typical energy loss time. As protons approach the reverse shock, they can interact with 
photons via photomeson production, if proton energy is above the threshold for photomeson production. It turns out 
that also in this region the photon density is high enough to ensure that almost all protons that have sufficiently high 
energies are depleted due to photomeson production. However, neutrino energy from this region is smaller than that 
from the pre-decelerated jet. 

The pp interaction also occurs in the shocked jet plasma, and a further neutrino source is due to non-thermal proton 
interactions with cold protons in the star after they escape the shocked jet plasma. For more detailed predictions, 
these processes should be taken into account. Yet, given the high efficiencies of the photohadronic processes for our 
typical parameters, we expect them to be at most subleading at the high energies we are mostly interested in. 

Finally, note that synchrotron, IC, and adiabatic cooling of pions and muons are i mportant in shaping our 
final result. The trea tment of these cooling processes is also similar to that used in (jMurase fc Nagataki"^ : 
iMurase fc Nagataki '06f) . Analytic considerations on effects of these cooling processes can be found in 
(jRazzaoue et al. '03f ). 
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